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Abstract. The aim of this paper was to describe changes in the spatial distribution of living trees in planted and 
untended birch stands, as well as to determine the spatial differentiation of the diameters and height of their constituent 
trees. Analysis was performed in 9- and 39-year old stands. 

The structural diversity of birch stands was described using the following spatial indices: Clark-Evans index 
(CE), spatial differentiation index for DBH ( TD ) and height (Til). The CE index describes the horizontal distribution 
of trees within a stand and the TD/TH indices indicate the size differentiation between neighbouring trees within 
a stand (DBH and H, respectively). The indices mentioned above concern small spatial scales. Our results indicate 
a regular spatial distribution of living trees on the ground at age nine, shown by high CE index. After 30 further 
years of stand development CE values were distinctly lower and no longer significantly different from the value for 
a random distribution of trees. At age 39, the 2 stands had a random spatial distribution and in the planted stand living 
trees still retained a regular spatial arrangement. Analysing the spatial differentiation of diameter and height of trees 
revealed that there was little differentiation between neighbours in forest stands in neither their diameter and height. 
The highest average values for the TD index decreased for stands during the 30 years since their initial measurement at 
age nine. A similar trend was observed for height differentiation (TH ), which decreased from a maximum at age nine. 

Key words: spatial structure of stands, size differentiation of trees, silver birch, plantations, semi-natural forest 



1. Introduction 

Due to the longevity of trees, a stand is subject to 
two types of variation: cyclical, relating principally 
to growth, and continuous, occurring over the very long 
period of a tree stand’s life and primarily related to 
developmental processes (Szymanski 2000). The most 
important of these include the natural differentiation of 
trees positions in the stand layers (biosocial position 
of trees), resulting from the processes of their natural 
development, the natural process of self-pruning 
and the process of stand ageing. To these should 
be added the creation of the stand-specific spatial 
structure, resulting from various factors, biotic as well 
as abiotic nature. Natural processes in a tree stand 
occur with varying intensity during specific stages of 
development, affecting its dynamic. Understanding 



the population dynamics of trees in artificially planted 
forests - whose share in the total area of forests both 
in Poland and in the world is significant - seems 
particularly important from the point of view of forest 
sustainable management. 

Until recently, tree stand structure was described 
most frequently with the use of such parameters 
as average size (thickness, height, etc.), average 
density or cross-sectional diameter at breast height. 
Currently, the spatial aspect is used to describe stand 
structure with increasing frequency, relating to the 
horizontal distribution of trees in the stand, spatial 
size differentiation or the spatial mixing of species in 
a multi species forest (Gadow, Hui 2002; Aquirre et 
al. 2003; Boyden et al. 2005; Szmyt, Korzeniewicz 
2010). Of particular importance are changes in 
stand structure occurring during its various stages of 
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development. The long life of trees, however, means 
that studying forest dynamics requires long-term work 
and analyses, repeated every several or several dozen 
years (Szymanski 1964; Ward et al. 1996; Kenkel et 
al. 1997; Hessburg et al. 2000; Woods 2000; Szmyt 
2004; Wolf 2005; Saunders, Wagner 2008; LeMay 
et al. 2009). Most research on forest dynamics has 
been conducted on natural forests; much less so on 
artificially established forests. The fact that the vast 
majority of forests in Poland and Europe originate 
from artificially planted stands is precisely the 
reason why research needs to be conducted on such 
populations. This will enable comparisons to be made 
of the structural variations between managed stands and 
natural forests, as well as to determine the influence of 
forest management on the development of the mentioned 
populations (Commarmot et al. 2005; Bilek et al. 2011). 

Managed forests are perceived as wood lands with a 
simplified structure, in the spatial aspect also however, 
the dynamics of these populations are relatively poorly 
understood. Supplementing this knowledge may lead to 
the development of new rational management principles 
that take into account naturally occurring ecological 
processes. Skilfully undertaken management activities 
may lead to an increased structural diversity of stands, 
which in the case of artificial populations is especially 
important, not only ecologically but also economically 
(Latham et al. 1998; Lahde et al. 1999; Lexerod, Eid 
2006). One of the factors determining the dynamics of a 
tree stand is the way the trees occur on the surface area, 
as this determines the influence of various ecological 
processes (e.g. competition, mortality). On the other 
hand, the same processes can shape the spatial structure 
of the stand in the next phases of development (Leps 
1990; Stoyan, Penttinen 2000; Comas, Mateu 2007, 
Eichhorn 2010). 

The present paper is on birch forests, initially planted 
in regular rows and left untended. For this reason, these 
are as provide good material to study the dynamics of 
artificially planted, untended tree stands. 

The purpose of this paper is to investigate the 
dynamics of the spatial structure of birch stands during 
the intensive growth stage of young trees, thus in a period 
when proportions are rapidly changing, in particular, 
the analysis included: 1) the placement of live trees and 
b) the spatial differences in tree height and thickness. 
The spatial variations of tree stands are defined with 
the help of indices of structural variation that can be 
used in forestry practice, such as in performing a forest 
inventory. 



2. The study subject and methodology 

The study was conducted at the experimental plots 
of the Department of Silviculture, Faculty of Forestry 
of the University of Life Sciences in Poznan, founded 
in the spring of 1970 by Prof. Stanislaw Szymanski. 
Administratively, the study area is located inthe 
southern part of the Wielkopolskie Voivodeship in 
K^pinsk County at the Forest Experimental Station of 
the University of Life Sciences in Poznan within the 
Wielislawice Forestry District in sections 25 hand 25i. 
The geographical coordinates of the area are: = 18° 

03 ’E, q)=51° 12’N. 

According to the regionalisation for natural and 
forested areas adopted in the country, the site is part of 
the Krainy Slqskiej (Land of Silesia), Wroclaw District, 
Rownina Olesnicka (Olesnicka Plain) mezzoregion 
(V.2.g) (Trampler et al. 1990). The average rain fall 
for the year is in the range of 550-600 mm (Schmuck 
1959; Ceitel, Wawro 1999). The birch stand selected for 
the study are growing on post-agricultural red podzolic 
soils formed from loose glacial sands, in all three areas 
(marked in the text as I, II and III), fully prepared soils 
and similarly aged planting stock (1-year old seedlings) 
were used with an initial spacing of 1. Ox 1.0 m. To date, 
the stands were left untended. 

Methodology 

In order to analyse the spatial differentiation of birch 
stands, the following measurements were taken during 
field work: root collar thickness of 10-year-old young 
tree stands (cm), diameter at breast height of older stands 
(D , cm) and total height (II. m) of all live trees. The 
position of each tree was described using rectangular 
coordinates. Measurements and observations were 
carried out in the autumn of 1978 and 2009 at three 
square areas with a side length of 20 m. The following 
spatial indices of structural diversity of the stands were 
calculated: the Clark-Evans index of aggregation (R), 
spatial differentiation index for diameter at breast height 
(TD) and spatial differentiation index of height (Til). 

The Rindicator enables the type of horizontal 
distribution to be specified of the studied population, 
based on the distance between the nearest neighbouring 
trees in a stand (Clark, Evans 1954; Brzeziecki 2002; 
Kint 2004). The method is based on comparing the 
average distance between nearest neighbours in the 
test stand and the expected (theoretical) average 
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distance between them at their random distribution. The 
mathematical formula for the R index with Donnelly’s 
modifications is as follows (Kint 2004; Szmyt, Ceitel 
2011 ): 

if r. 

r JvZj ' 

R=-L= fid 

r E 0,5(|) 1/2 + 0,0514^+0,041^ 

where: 

r 4 - average distance between nearest neighbours, 

r — average expected (theoretical) distance between 
nearest neighbours, 

r— distance between an individual tree and its nearest 
neighbour (m), 

N - number of individual trees in the measured area, 

A - size of the area (m 2 ), 

P - diameter of the area (m). 

Tree stand with a random distribution of trees 
is characterised by an index of R=l; when trees are 
arranged in clusters, then R< 1. The value of R> 1 
indicates an even distribution of trees in the test stand. 
The significance of deviations of the calculated index 
from the theoretical pattern can be determined using the 
standard Z normal distribution (Kint 2004): 



The index value is calculated for each tree in the 
stand, and then averaged to indicate the average spatial 
differentiation of the size of the trees within the stand. 
The index can have values from 0 to 1, with 0 indicating 
no spatial variation in thickness and that the trees are all 
the same size, where as the value of TD—l indicates a 
very high variation in the dimensions of the trees (Kint 
2004; Bilek et al. 2011; Szmyt, Ceitel 2011). 

In the case of the spatial differentiation index values 
( TD and Til), the following class differentiations were 
adopted: 0.00-0.20, very small variation (nearest 
neighbour difference in size is not more than 20%); 
0.20-0.40, small variation (difference in size is 20- 
40%); 0.40-0.60, medium variation (difference of 
40-60%); 0.60-0.80, high variation (difference of 
60-80%); and 0.80-1.00, very high thickness variation 
(difference of 80%) (Aguirre et al. 2003; Kint 2004; 
Vorcak et al. 2008). 

All calculations in this study were carried out using 
the SIAFOR ver. 1.0 software (Kint 2004). 

3. Results 

The horizontal distribution of the trees 



c = - 



a 



r E 



a 



r E 



0,26136 



where: 

a - the standard deviation/- fora random distribution 

rE E 

of density p. 



The index of spatial diameter and height 
differentiation of the trees (TD and TH, respectively) 
determines the similarity in the diameter or height of 
trees in the stand at the smallest spatial scale, i.e. the 
scale of the immediate neighbour. The mathematical 
formula for the index involving three nearest neighbours 
is denoted by (Pommerening 2002; Aguirre et al. 2003; 
Kint 2004; Szmyt, Ceitel 2011): 

min ( s ,Sj ) 
max (5, ,Sj ) 




where: 

i - indicator for the reference (central) tree, 
si - size (DBH or height) of tree i, 
sj - size (DBF! or height) of tree j. 



The value of the aggregation index R indicates a 
significant change in the horizontal distribution of live 
trees over a 30-year period in all three stands. The initial 
value of the index (for trees at 9 years of age) clearly 
indicated a very regular distribution of individuals in 
each stand (R> 1, =0.05). Specific stands of this age 
were characterised by a similar value of the Clark-Evans 
index (Table 1). During 30 years of further development, 
the index value of these stands decreased significantly. 
In their 39 th year of life, a random distribution of live 
trees was determined for two stands (R=1.05 and 
R= 1 .00), while trees in another stand continued to be 
regularly distributed (R=1.19). 



Table 1. Average value of spatial diversity indices of 
distribution ( R ), DBH differentiation (TD) and height 
differentiation (TH) in untended birch stands at different ages 



Forest 


9 years 


39 years 


stand 


R 


TD 


TH 


R 


TD 


TH 


I 


1 75 ** 


0.30 


0.22 


1.05 


0.26 


0.23 


II 


1.72** 


0.37 


0.26 


\ 19** 


0.29 


0.17 


III 


1.73** 


0.39 


0.30 


1.00 


0.29 


0.20 



** level of significance a=0.0 1 
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Table 2. Statistical characteristics of DBH (D) and height ( H) of live trees in untended birch stands at 9 and 39 years of age 



Forest stand 


Tree 

age 


Number of 
trees 


G„ 

ha 


Survival 


Average 


Min 


Max 


Standard 

deviation 


Coefficient 
of variation 


years 


m 2 


% 


D 

(cm) 


H 

(m) 


D 

(cm) 


H 

(m) 


D 

(cm) 


H 

(m) 


D 

(cm) 


H 

(m) 


D 

(%) 


H 

(%) 


I 


9 


323 


20.6 


100.0 


5.4** 


4.2 


0.5 


0.8 


11.3 


6.4 


1.84 


10.92 


33.54 


25.79 


II 




312 


23.4 


96.6 


5 7** 


4.7 


0.6 


2.2 


12.1 


8.3 


2.46 


14.42 


43.41 


30.51 


III 




314 


17.2 


97.2 


4 g** 


4.1 


0.6 


4.0 


10.8 


7.6 


2.23 


14.55 


46.73 


35.20 


I 


39 


48 


13.7 


14.9 


11.8 


17.1 


7.0 


4.5 


18.0 


21.8 


2.71 


3.92 


23.06 


22.85 


II 




45 


17.9 


13.9 


13.7 


20.2 


8.0 


14.4 


21.0 


25.7 


3.73 


3.10 


27.23 


15.32 


III 




38 


16.4 


11.8 


14.1 


20.5 


7.0 


12.8 


24.5 


28.0 


4.33 


3.84 


30.55 


18.73 



- root collar thickness 



% % 




Figure 1. Distribution of the spatial differentiation index for DBH (TD) in particular classes of differentiation at different ages 
in untended birch stands 



Spatial variation of tree diameter 

The average rate of variation in the diameter of 
trees in the stand exhibited little diversity from the 
very beginning, with trees not differing more than 40% 
from their nearest neighbours (Table 1). Over 30 years, 
this difference significantly declined even further in all 
three stands. The reduced variation in thickness was 
also confirmed by the coefficient of variation of this 
characteristic, whose values also significantly decreased 
in the 39-year-old stands (Table 2). 

Figure 1 shows the proportion of trees in the adopted 
classes of diameter differentiation. The largest number 
of trees in the studied stands at 9 and 39 years of age 
were in the class of 20-40% diversity. At the beginning 
of the study period, a distinct proportion of trees were 
in the medium class of 40-60% diversity in two stands, 
whereas at the end of the study period, the proportion 
of trees in this class decreased in favour of the class 



with the least diversity of diameter between immediate 
neighbours. The figure shows that in the 9 th year of 
age, only the thickness of a small number (%) of trees 
differed significantly from their nearest neighbours 
(difference of 60%). After the stands developed for 
a further 30 years, trees exhibiting such a difference in 
diameter at the breast height were not found. 

Spatial variation of tree height 

Table 1 shows the average value of the spatial 
diversity index of tree height for individual stands. 
It shows that tree height variation in both the 9- and 
39-year-old stands is small and decreases with age. Only 
in Stand I the average index value increased slightly 
from 0.22 to 0.23 at the end of the analysed period. 
The decrease in height was also confirmed by the 
coefficient of variation calculated for the height of the 
trees (Table 2). 
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% 




Figure 2. Distribution of the spatial differentiation index for height (TH) in particular classes of variations at different ages in 
untended birch stands 



The proportion of trees in each class of height 
differentiation in comparison with its neighbours is 
shown in Figure 2. Both at the beginning and at the end of 
the analysed period, the highest number of trees belonged 
to the class of up to 40% of height diversity. With age, 
the proportion of trees in the class with the lowest 
diversity increased in two stands, at the expense of the 
other classes. The figure also shows that the proportion of 
trees significantly differing in height compared with their 
neighbours decreases with age, reflecting a reduction in 
the height variation of the tree stands. At the age of 9, 
the proportion of young trees differing in height by more 
than 40% ranged from 8.9% to almost 22%. At the end of 
the study period, only one stand exhibited this variation 
(40%) in about 8% of its trees. 

4. Discussion 

An initial regular spatial distribution of trees 
is typical of untended stands planted in a specified 
pattern, and the results confirm previous studies. The 
regular spatial distribution of live trees in plantations 
and young stands is confirmed by, for example, Szmyt 
(2004), who analysed the spatial structure of artificially 
planted, untended pine, spruce and oak established in 
various patterns of rows. Regardless of the pattern used 
in planting, all the analysed stands exhibited a distinctly 
regular spatial distribution at the time of harvesting, 
which resulted from the initially imposed pattern. 
Flowever, this need not be a permanent feature, as it 
may vary in subsequent phases of stand development. 
The direction of these changes in artificially planted 
stands is not sufficiently understood and requires further 
studies (Szwagrzyk 1992; Sekretenko, Gavrikow 1998; 
Szmyt 2004; Szmyt 2010). 

In the case of naturally regenerated forests, the 
dynamics of tree distribution is much better known 



and most often starts with a clustered distribution 
in the regeneration phase, then goes to a random 
distribution and finally a regular distribution in older 
age-class forests (Gavrikow, Stoyan 1995). In the early 
developmental phase of artificially planted tree stands, 
the clustering of individuals should not be expected, and 
this was confirmed by our research. However, in older 
stands, when the effect of the initial spacing lessens due 
to natural processes, the future direction of changes can 
occur as it does in natural forests (Szmyt 2004; Szmyt, 
Korzeniewicz 2010). A random distribution of live 
trees in older stands of artificial origin was observed 
by Szwagrzyk (1992) and Szmyt (2010). However, 
Sekretenko, Gavrikow (1998), Pommerening (2002), 
Kint et al.(2003), Mason et al. (2007) and Szmyt, Ceitel 
(2011) found a regular distribution of trees in stands of 
similar origin. 

The reasons for the regular distribution of trees in 
stands are seen primarily as the result of competition 
between individuals, but competition is not always the 
most important process occurring in populations, as often 
many other factors and ecological processes play a more 
significant role in the formation of the characteristic 
spatial structure (Kenkel 1988; Szwagrzyk 1992; Murrell 
2009). The impact of competition resulting from the 
high density of individuals in naturally regenerating tree 
stands is especially visible in natural forests, when the 
intensity of individual growth is greater in areas with a 
higher density, often leading to the random or regular 
distribution of the remaining trees. In artificially planted, 
managed stands, the choice of initial planting pattern and 
density may limit competition among individuals. As 
Szmyt (2004) noted, the sparser the initial plantings, the 
more durable the regular distribution of trees in the stand. 

Often, single species and similarly aged stands are 
considered to have little structural differences. Our 
results support this view in the diversity of both the 
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thickness and height of the trees. In untended birch 
stands, the coefficient of variation of thickness and 
height, as well as the spatial diversity indices of both 
these characteristics, indicate little variation throughout 
the 39 years of the analysed period. This diversity was 
greater in stands of young birch trees and significantly 
decreased with age. During cultivation, before the tree 
canopy closed, the growth of young trees was individual 
in nature and residted from the interaction of genotype 
and environment. As soon as the tree canopy closed, 
competitive interaction between individuals took on 
greater significance. As a result of this competition, the 
weakest individuals are eliminated, usually those with a 
small diameter at breast height and low height. 

The study results presented of the birch stands are 
consistent with other research on artificially planted 
stands (Pommerening 2002; Szmyt, Korzeniewicz 2010; 
Szmyt, Ceitel 2011). In examining the diameter variation 
of artificial pine stands of different ages, Szmyt (2010) 
found that the least variation occurred in the oldest 
stands (90-year-old), but it was also small in 43- and 
65-year-old stands. Little variation of diameters was 
found in the artificially planted pine stands analysed by 
Mason et al. (2007). Dro|3er and Liipke (2004) observed 
a similar trend in single species managed stands of 
beech aged 148 and 166 years. Younger stands (80 and 
120 years) were more varied in this feature. Little height 
variation was found by Kint et al. (2003) in artificially 
planted pine tree stands. As Brzeziecki (2005) noted, 
the diversity of tree stands depends not only on various 
ecological processes, but also on how they are managed 
(e.g. thinned). The author found that in a sample pine 
stand, diameter variation was small before thinning, 
while thinning from below, thus preceding the natural 
removal of trees, decreased variability even more. On 
the other hand, selective thinning, the recommended 
practice, promotes the formation of greater spatial 
variation in a stand. 

Crecente-Campo et al. (2009) analysed the impact of 
thinning on the structural diversity of 50- to 60-year-old 
artificially regenerated pine stands. The control plot (no 
thinning) had a random distribution of trees, and indicators 
of the structural diversity of thickness and height decreased 
as a result of the natural elimination of trees from the 
stand. The influence of thinning on the aforementioned 
characteristics depends, among other things, on its intensity 
(Crecente-Campo et al. 2009). Similar conclusions were 
reached earlier by Pretzsch (1999) who analysed the effect 
of different management practices on the diversity of 
mixed spruce and beech stands. 



5. Conclusions 

The results obtained allow the following conclusions 
to be proposed. 

Despite the initial pattern of distributing seedlings 
during planting, horizontal distribution can become 
random as a result of the natural processes eliminating 
trees from the stand. If left untended, this processis the 
main one influencing the horizontal distribution of trees 
in a stand. 

The results presented confirm the thesis that the lack 
of tending in birch stands during the initial stage of their 
development leads to a decrease in thickness and height 
diversity. Based on the natural ecological processes 
occurring in birch stands, increasing their diversity is 
not possible. 
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